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This paper describes a simple and general method for deriving the inelastic collision term in 
the electron Boltzmann equation for scattering from a coupled electron-phonon system, and applies 
the method to the case of doped polar semiconductors. In the Born approximation, the inelastic 
differential scattering rate VK'"*' can be expressed in terms of the nonequilibrium total dynamic di- 
electric function, which includes both electronic and lattice contributions. Within the random-phase 
approximation W"^"^' separates into two components: an electron-electron interaction containing the 
lO I nonequilibrium distribution function for excitations of the electron gas, and a Frohlich interaction 

0^ ■ including the phonon distribution function and self-energy due to polarization of the electrons. Each 

of these two interactions is screened by only the electronic part of the total dielectric function which 
contains the high frequency dielectric constant, unlike commonly used expressions which contain 
■ the static dielectric constant. Detailed balance between plasmons and electron-hole pairs in steady 

state is used to eliminate the nonequilibrium plasmon distribution from the Boltzmann equation, 
resulting in a dynamically screened electron-electron collision term. The phonon self-energy modi- 
0^ ' fies the longitudinal optical phonon dispersion so that two hybrid normal modes contribute to the 

electron-phonon collision term. 
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. The electron Boltzmann equation for inelastic scattering in solids is usually established either by making approxi- 
mations to more general equations derived with powerful nonequilibrium Green's function techniques or simply 
by ansatz with the aid of analogy to simpler physical systems. An alternative method described here proceeds by 
using time-dependent perturbation theory to determine the inelastic collision term from the total dynamic dielectric 
function for the nonequilibrium coupled electron-phonon system. This paper treats the case of a system of polar 
optical phonons in the presence of conduction electrons and Coulomb interactions, though the method is more gen- 
] erally applicable. Except for the nonequilibrium character of the Boltzmann problem, the approach is similar to Van 
I I Hove's correlation function method applied to doped polar semiconductors j6| to determine the inelastic lifetime 
'"O of a quasiparticle excitation of an equilibrium system. In addition to providing physical insight, the total dielectric 
, function approach is a systematic way to include scattering against coupled electron-phonon modes in the Boltzmann 
^ ' equation, so that the influence of mode coupling on mobility can be studied. Numerical results for mobility in n-type 
GaAs are presented in reference 0. 

In the Born approximation, the interaction between a conduction electron and a coupled electron-phonon system 
can be represented as an effective electron-electron interaction screened by the total dynamic longitudinal dielectric 
function, er(q, w), which includes contributions from both electrons and phonons P-JlO|. In this context, er describes 
■ the linear response of the electron-phonon gas not to an externally applied potential, but to the internal potential due 
to the probe electron plus the induced density fluctuations of the gas itself. The effective interaction has the form of 
a total density autocorrelation function. Unlike the lifetime problem the entire electron-phonon gas is away from 
equilibrium in the dc transport case, and the correlation function must be evaluated for the nonequilibrium ensemble. 
It follows that the effective electron-electron interaction is screened by the nonequilibrium total dielectric function in 
this case. In fact, taking account of the nonequilibrium character of the interaction is essential to obtaining the correct 
form of the inelastic collision term, which contains nonequilibrium distribution functions for the system's excitations. 

The effective electron-electron interaction can be separated into a purely electron-electron term and a screened 
electron-phonon interaction. In the random-phase approximation (RPA), the electron-phonon part contains a phonon 
self-energy which arises from the polarization of the electron gas [p^ . The self-energy correction modifies the lon- 
gitudinal optical (LO) phonon dispersion in doped polar semiconductors, producing two hybrid normal modes with 
phonon strength in both. This paper determines the inelastic differential scattering rate W™^^ for the nonequilibrium 
case in the RPA Born approximation. The result consists of an electron-electron component containing the nonequi- 
librium distribution function for excitations of the electron gas, and an electron-phonon component containing the 
phonon self-energy and distribution function. Each of these two interactions is screened by only the electronic part 
of the total dielectric function, which contains the high-frequency dielectric constant eoo- Thus, the RPA Born result 
differs from commonly used expressions for screened interactions in doped polar semiconductors which contain the 
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static dielectric constant eg. Examples using eg rather than eoo for n-type GaAs include a calculation of dynamically 
screened electron-electron scattering , as well as static screening approximations for electron-LO phonon scatter- 
ing [|l2| and electron-electron scattering For more strongly polar doped semiconductors (where eo 3> £00), it is 
essential to treat this aspect of the screening problem correctly. 

The collective modes have a finite lifetime for decay into electron-hole pairs. This Landau damping process provides 
an indirect mechanism for single-particle electron-electron scattering. For the electronic system in steady state, 
Landau damping creates a detailed balance relation between plasmons and electron-hole pairs which can be used to 
eliminate the nonequilibrium plasmon distribution from the Boltzmann equation in favor of the single-particle electron 
distribution. In the approximation of no collision damping, the result is an electron-electron collision term screened in 
the dynamic RPA, containing the effects of electron-plasmon scattering. For direct gap doped semiconductors where 
Umklapp processes are negligible, this collision term conserves total electron momentum and cannot by itself degrade 
an electrical current if energy surfaces are spherical. However, by rearranging the electron momentum distribution, 
electron-plasmon and electron-electron scattering can have a significant indirect effect on the outcomes of other 
scattering processes. 

The low-energy hybrid mode arising from the plasmon-phonon coupling in doped polar semiconductors gives in- 
creased low-energy electron-phonon scattering strength compared to the uncoupled case. Numerical calculations 
including the plasmon-phonon coupling have given shorter inelastic lifetimes |^,^^ and enhanced hot electron energy 
relaxation |15 1^, especially at low t emp eratures and low doping. A simple way to treat the coupled modes is to use 
the plasmon-pole approximation |p7| , |l8| on the phonon self-energy. This procedure is used here to include scatter- 
ing against the hybrid modes in the electron-phonon collision term so that mode coupling can be treated with the 
Boltzmann equation in a computationally efficient manner. 

This paper is organized as follows. In section II, W™'''' is derived within the Born approximation in terms of the total 
nonequilibrium dielectric function for the coupled electron-phonon system. The RPA is invoked to separate VJ^™"^' 
into electron-electron and electron-phonon components. Section III deals with the electron Boltzmann equation by 
eliminating the nonequilibrium plasmon distribution from the electron-electron collision term. The electron-phonon 
collision term is simplified by applying the plasmon-pole approximation to the phonon self-energy. 



II. INELASTIC DIFFERENTIAL SCATTERING RATE 



A. Electron-phonon gas model 

Consider a probe electron scattering from an interacting gas of electrons and polar optical phonons. Interactions 
with impurities and acoustic phonons can be added by hand later. The probe electron is treated as a quasiparticle 
with states |k), effective mass m*, and dispersion relation Ek — h^k^ /2m* . The Boltzmann equation relies on the 
validity of the quasiparticle description by assuming that the energy-momentum relation remains well defined when 
collisions are included. Interactions with the electron-phonon gas produce transitions between the quasiparticle states, 
but do not broaden them significantly. 

For direct-gap polar semiconductors, local field effects are unimportant so that electrons couple only to the macro- 
scopic electric field set up by the LO phonons, and do not interact with the transverse optical (TO) phonons. The 
Hamiltonian for the electron-LO phonon gas is 

Hg - 5I^fe4,.Ck,. + ^?iwTo(<aq + J) + ^IZ^r \pU^)pt{(1) ~ n] , (1) 
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where 



vr = ^ = ^, (2) 

£00 Q £00 

and where the electron and LO phonon annihilation operators are denoted by Ck and Oq, respectively, a labels the 
electron spin, and ujto is the TO phonon frequency. The relevant phonon frequency in Hg is the bare one, ujto- The 
ion- ion Coulomb interactions explicitly included in Hg shift ujto to the longitudinal optical frequency wloi so that 
only lolo appears in expressions for observables |lO| . The first and second sums in Hg describe the noninteracting 
systems of electrons and phonons, respectively. The last sum includes electron-electron, electron-phonon, and ion-ion 
Coulomb interactions. The total density operator pT is the summation of the density operators of the electrons and 
LO phonons. 
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M,) = n-./Mc)H-iJ,(^)" 



Ml), (3) 

where A{q) — aq + oL^, M is the reduced mass of the basis ions, Ze is the effective ionic charge, n is the electron 
density, f2 is the system volume, and the Fourier components of the electron density operator are 

p(q)=^e-''i'-'=5:ct+q,,Ck,.. (4) 

i k,(T 

The electron self-interaction has been removed from Hg by subtracting the term with N ^ the number of electrons in the 
system. The influence of interband electronic transitions is taken into account by screening the Coulomb interaction 
Vq with the high-frequency dielectric constant too- 

The effects of screening by conduction electrons and LO phonons are contained in the interaction term in Hg and 
will appear explicitly in the total dielectric function below. The RPA or mean-field approximation for the interaction 
term is made by working with eigenstates |7i) of the Hamiltonian Hg which neglects the terms in Hg that are quadratic 
in density fluctuations pT — (pt) 

Hg\n) = En\n) (5) 

Ha = Hg- 1 ^«-[4(q) - (pt,(q))] [pj.(q) - (p^(q))] . (6) 
q 

Hg includes Coulomb interactions at the Hartree level. Excitations of the electron-phonon gas described by Hg include 

both electron- hole pair and collective excitations. 

Define P{n, k ^ m, k — q) as the probability per unit time for the probe electron to make a transition from |k) to 
k — q) while the electron-phonon gas makes a transition from \n) to |m), assuming that initially |k) is occupied and 
k — q) is unoccupied. The interactions included in Hg produce charge density correlations, resulting in an interaction 

potential between gas excitations and probe electron which is screened and assumed to be weak. Using Fermi's Golden 

rule (or the Born approximation) with the probe particle plane wave states |k), the transition probability is 

P(n,k~>m,k~q) = — f-^j \ {m \ pt^iq) \ n) \^ S{E„ - + nu;k,k-g) (7) 

= [-^) J rfte-'=.-'*(n|pT(q,t) l™)(m|p^(q,0) |n) (8) 

PT{q,t)^e^"^'/^pr{q)e-^"^'/\ (9) 

where LOk,k-q = - (k - q)^]/2m*. 

The inelastic differential scattering rate VF™*^' (k, k — q) is the probability per unit time for an electronic transition 
from the occupied state |k) to the unoccupied state |k — q) due to scattering from the electron-phonon gas. It is 
found by summing P{n, k ^ to, k — q) over all possible final states | to) of the gas and averaging over its initial states 
n). For the mobility problem, the average must be taken over the nonequilibrium gas ensemble. 



B. Nonequilibrium ensemble averages 

Consider that an electric potential f/(r, t) is applied to the system of probe electron plus electron-phonon gas. This 
influence is included by adding the term 

He,t{t) ^ejd\ U{r, t)pT{r, t) (10) 

to the Hamiltonian of the system. H^xt is independent of time in steady state, but it is useful to require the potential 
to be switched on at a time in the distant past when the gas was in thermal equilibrium. This requirement makes it 
possible to use the Kadanoff and Baym method p9| ] for evaluating a nonequilibrium expectation value. This method 
uses the equilibrium density matrix to evaluate an ensemble average at a time after U has been switched on by using 
the interaction picture to include the dependence on H^xt explicitly. The nonequilibrium expectation value of an 
operator 0(r, t) is ^ 
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{0{r,t))„ = {Ou{r,t)) = Z~^Tt {e-^''^Ou{r,t)} (11) 
Ou{r,t) = V-'(t)0{r,t)V{t), (12) 



V(i) = rexp 



—I 



(13) 



where ( ) denotes an average over the equilibrium population of the electron-phonon gas, /3 — l/fesT, Z = Trje"'^^"} 
is the partition function for the equiUbrium gas, and T denotes the time ordering operator which orders operators 
with earhest times to the right. 

In the presence of t/(r, t), the transition rate P(ri, k — » m, k — q) in the Born approximation is unchanged from its 
equihbrium form except that the influence of He^t must be included in the time dependence of the matrix element. 
Using equations (ll]-13|) to determine W*"'^'(k, k — q) from P(n, k ^ m, k — q) for the noncquilibrium case yields 

T^-^'(k,k-q)= f^j NSTu{c\,UJk,k-q), (14) 
1 1"^ 

STu{ti,uj)^- dte'"*(pTC/(q,t)PTO(-q,0)) (15) 

J —OO 

PTuici, t) = V-\t) e'"^'/'' pT e-'^^*/" V{t), (16) 

where Srui^T^) is the noncquilibrium total dynamic structure factor or spectral function for pruj describing the 
density fluctuation excitations of the nonequilibrium electron-phonon gas. Similarly, the time-reversed rate is 

W^"^'(k-q,k)= f^j NSTu{-ci,-iOk,k-,). (17) 

When the gas excitations are not too highly damped so that they may usefully be regarded as elementary excitations, 
and when the deviations from equilibrium are slowly varying in space and time on a scale determined by the frequencies 
and decay rates of the elementary excitations, it is possible to define nonequilibrium distribution functions and response 
functions 0J^. The imaginary (dissipative) part of the response function for the nonequilibrium coupled electron- 
phonon system is 

-1 

Im[xTt/(q, ^)]^^J dt e'-*([pTt/(q, t), PT[/(-q, 0)]) (18) 

— 77 

= — [Stu{ci.uj)~Stu{-ci.-oj)] (19) 

where [p, p'] is the commutator pp' — p' p. 

In general, there is no simple relation between a nonequilibrium spectral function and corresponding response 
function like the fluctuation-dissipation theorem (FDT) for the equilibrium case: 

5T(q, UJ)^ — {N\uj) + 1] Im[xT(q, ^)\ , (20) 
n 

where N^(lo) = (^f^^^ _ 1)^^ is the equilibrium Bose-Einstcin distribution function, while St and Im^T are the 
equilibrium forms of ( p^ and (|l8|). In particular, since the relation between the time- reversed nonequilibrium spectral 
functions Stu{(Ij^) ^ind STui~^i is not known, Im[xT[/(q, t^)] is not determined by S'T[/(q, w) alone. It is 
possible, of course, to define unknown dimensionless functions iV>(q, w) and A^<(q, lj) such that 

5TC/(q,w) = ^^>(q,w)Im[xTO(q,w)] (21) 
STui-ci,-io) - ^iV<(q,c^)Im[xTC/(q,t^)]. (22) 
Equations (^and (|l|)-(||) then imply 

N>{q,Lo)~N<{q,Lo) = l. (23) 
Adopting the notation N^{q,uj) = N{q_,Lu), we have from (p^),(21 ), and (p3|). 
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W™^'(k,k - q) = [N{q,u;k,k-q) + 1] lm[xTu{<l,0Jk,k-q)] ■ (24) 

Equivalently, the Born aproximation for W™'^^ can be expressed in terms of the imaginary part of the total screened 
Coulomb interaction Vq/exu between the probe particle and the electron-phonon gas, 

H^™^'(k,k - q) - ^ [N{ci,UkM-,) + 1] Im [e^]j{ci,UkM-,)\ ■ (25) 

This follows from the relations between the total susceptibility XTu{^.^'^)^ polarization P7'(7(q, w), and dielectric 
function eTu{''\i^) for the nonequilibrium coupled system, 

I X fT(7(q,w) 

XT(7(q,w =^1 —B — 7 7 (26) 

1 -^^^^(q,^) 

eTC/(q, ^) = eoo [1 - v^Ptu[ci, lo)] . (27) 
The relation (u^)^ImxT(7(q, w) — VqluYe^]j follows from ( p6| ) and (p7|), and yields ( p5| ) from (| 



The weight factor iV(q, w) in (25) plays the role of the distribution function for the nonequilibrium excitations of 
the coupled electron-phonon system. Just as for the equilibrium case in equation ([20|), the structure factor has been 
written as a product of a function (Im[ey^]) describing the strength of interactions with the electron phonon gas, and 
a function giving the occupation probabilities for the gas excitations. However, it should be emphasized that the FDT 
has not been used to derive (|25| ) and N[(\,uj) is an unknown function. 

C. Separation of the total interaction 



The total nonequilibrium structure factor Stu defined by (jl^) is exactly separable into a purely electronic part Sfj 
plus the remainder S^^ which includes electron-phonon and ion-ion Coulomb interactions. Using the nonequilibrium 
version of the total density operator (||), 

5TO(q,c^) =5&(q,c^) + ^&'(q,^), (28) 
5&(q,c.) = — y dte"^*(pt;(q,t)pt;(-q,0)), (29) 

Sfj\q,io)^- I dte^-'{Zq[ ^J ^ ] (pc/(q,t)Ac/(-q,0)) (30) 



N 1 ^\2MujTO^J 

+ (At;(q,i)pt;(-q,0))] + (Zg)2^^^!^(^t/(q,t)Aa(-q,0))| ■ (31) 

This section relates and S'p'' to corresponding electron and phonon parts of the response function, just as (|2l| ) 
relates the total structure factor to the total response function. To do so, equation (^5|) must be split into electron- 
electron and electron-phonon parts. Imposing the RPA for the dielectric function is a crucial step towards this goal. 
The response of the coupled system in the RPA is the sum of the electron and ion response taken separately |^,^ . 
In this approximation, the total screened interaction Vq/eru is the sum of a purely electron-electron part and an 
electron-phonon interaction that includes a phonon self-energy due to the polarization of the conduction electrons 

ill- 

ETC/ I £(7 r 

The electron-electron part of the total interaction is u^/ejy, where eu is the nonequilibrium RPA dielectric function 
determined by electron-electron interactions only. The polarization for the noninteracting electron gas P^^^ has the 
same Lindhard form as the equilibrium case, except that the nonequilibrium electron distribution function / is used 
instead of the Fermi-Dirac function |^ . 

eu{ci,u) = l-v^P^'\ci,u) (33) 
\l niOp^p+q ~ nu) — 10 
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The remaining part of the total interaction is a product of the screened electron-phonon matrix element M^|e£/| ^, 
and the nonequilibrium phonon Green's function Dif{q,u!) containing the self-energy correction M'^D^^^ xu to the 



noninteracting Green's function D 



(0) 



D 



(0) 

u 



Du[ci, Lo) = ; . (36) 



l-M2Z?^'^Xc/(q,w) 

> > 

Now define unknown dimcnsionlcss functions and N^f^ to relate the electron and phonon components of 
^^[/(q, Lj) in ( |2^ ) to the respective components of the imaginary part of Vqe^^ in (p2[). 



{v?^rSh{±<l, ±^) = v^N<ici,Lu)Im[e^\<i,u;)] (37) 



2 



{v^)'S^J'{±<l,±u;) ^ ^^7V<(q,c.)Im[i?^(q,^)]. (38) 
n I cjj I 



Finally, using (0) , (|2|) , (|3) , and (| 



iy™^'(k,k- q) = ^\v^[N,ici,Lu,^k-,) + 1] Im [e^'{q,u;k,k-,)] 



Af2 ] 

+ I f \l2 Wph{<l, ^k,k-q) + 1] Im [^[/(q, UJk,k~q)] ) ■ (39) 

I eu[q,i^k,k^q) \ 



Within the RPA, the expressions (H) and (Ig) for W^*"^' are equivalent. The last expression contains two unknown 
functions, Ng and Nph, playing the roles of the nonequilibrium occupation functions for the excitations of the electronic 
and lattice components of the coupled system. Notice that must describe a distribution with boson properties. 
This is because the density fluctuation produces excitations of the electron gas which conserve the total number 
of electrons, such as electron-hole pairs and plasmons, which have boson properties. 



III. ELECTRON BOLTZMANN EQUATION 

The time rate of change of the nonequilibrium electron distribution /(k) due to collisions is determined by the 
differential scattering rate M^(k, k — q); i.e., the probability per unit time for the transition |k) — ^ |k — q). The 
depletion rate of /(k) is just W(k, k — q) weighted by the probabilities that |k) is occupied and |k — q) is unoccupied, 

summed over all possible final states. The collision term, — {/ (k)}coz;, in the Boltzmann equation is this depletion 
rate minus the replenishment rate determined similarly. The electric field in the sample is F(r, t) = ~V[U{r,t) + 
{Uind{T^i i))]) where (Uind) is the average induced potential resulting from the system's response to the applied potential 
U. Thus, for a homogeneous system under the influence of a static electric field F, 



-eF 9/(k) 



h 9k L J coll 



i Y.{w{k, k - q)/(k)[l - /(k - q)] - W{k - q, k)/(k - q)[l - /(k)]}. (40) 



W^k, k — q) is a sum of rates for elastic and inelastic processes, W — W^^ -\-W™'^^ . Using the separation of the inelastic 
rate Wmei in (p9|), the inelastic contribution to the collision term separates into electron-electron and electron-phonon 
collision terms, 

{/}"?/ = • (41) 
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A. Electron-electron collision term 



The purely electronic component of W™'^'' in (^9|) gives the collision term 

- {/ {^)Yclu = 

-^E^rj [NMM-,) + l]Im[e^i(q,c.fe,fc„,)]/(k) [1 - /(k - q)] 
q I 

-7Ve(q,c^fe,fc_,)Im[e^i(q,L^fe,fc_,)]/(k - q) [1 - /(k)] | . (42) 

The next step towards putting the collision term into a useful form is to eliminate the unknown function N^, which gives 
the occupation probabilities of the density fluctuation excitations of the electron system. This step is straightforward 
if an association is made between N^. and the nonequilibrum plasmon distribution, denoted here by Np. In steady 
state, there exists a local detailed balance (LDB) between plasmons and electron-hole pairs which can be used to 
eliminate Nf. in favor of / in the electron-electron collision term |^,^. The LDB condition is equivalent to the 
FDT for the nonequilibrium ensemble [ p3| which applies when x determines the response not to external forces, but 
to internal fluctuation forces. It holds when the Landau damping rate 7 = 2ImP'^^^ is large compared to the time 
variation of /, so that the plasmons come into equilibrium with the electrons before f changes appreciably. For the 
mobility problem, / is stationary and the condition always holds exactly, except at small q values (where 7 vanishes) 
which do not contribute to the electron collision integral because of energy and momentum conservation restrictions. 
If LDB is used in combination with the "coUisionless damping" approximation equation (42) is equivalent to 



the standard single-particle electron-electron collision term, but screened with the dynamic RPA dielectric function. 
This may be shown easily for the case of weak damping, in which the plasmons are well-defined excitations, by 
considering the Boltzmann equation for the nonequilibrium plasmon distribution Np{q) = Np(q,ujq) with LUq = ^p{q) 
[^,|^. Plasmons are represented by the peaks in the spectral function Ime^^. In the RPA, this function contains no 
corrections for collision broadening, so that the relaxation of Np{q) is completely described by coherent absorption 
and emission of electron-hole pairs. 

Np (q) = ^ eI + ll/C^ + [1 - /(l^)] 



k 



-Npiq)f(k) [1 - /(k + q)] }Sinu;k,k+q - hujq) . (43) 



In steady state, Np= 0. Therefore, the relation between the distributions Np and / is necessarily 



Np{q)lm[P^'\ci,^q)] = —Sl'\ci,cOq) (44) 



lm[P^'\<l,u;)] = ^ E [/(P) - + '5(?i^p,p+, - hu;) (45) 
,(1)/ ^ 47r;i 



p 



p 



relating the response function Im[P^^''] and structure factor 5'^^'' for the noninteracting nonequilibrium electron gas. 

The detailed balance relation between plasmons and electron-hole pairss simplifies the electron-electron collision 
term (p2|), if we assume N^ — Np. Using the RPA interaction Im[e^^] — t;^|e[/|~^Im[P^^-'] and (^-46) gives 



if (kl)}coH = E I ^k-E,,,2 ^ + ^2 - Es E,)6{k^ + k2 - k3 - k4) 



{/(kl)/(k2)[l-/(k3)][l-/(k4)] 

-[l-/(kl)][l-/(k2)]/(k3)/(k4)}, (47) 
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when the notational changes k — > ki,p ^ k2 are made, the summation variable q is changed to ka = ki — q, and 
a summation over a delta function for momentum conservation is included. Notice that the nonequilibrium inverse 
dielectric function was needed to obtain the collision term with nonequilibrium distributions for all four (incoming 
and outgoing) electrons. The remaining factor of |e[/P in equation (|4^ ) appears only in its equilibrium form |ep if the 
Boltzmann equation is linearized with respect to the field strength F, since the other factors vanish in equilibrium. 
Also notice that eoo rather than eq appears in the electronic dielectric function eu{q,^) (see equation (|33|)) which 
screens the interaction in the electron-electron collision term (^), unlike some previous calculations using eo for 
electron-electron scattering in GaAs • 

It should be emphasized that (^^ includes the effect of electron-plasmon scattering through the resonance in e{q, lo) 
at Lo = ujq. For dc transport, the two-step process: 

electron — hole pair — > plasmon — > electron — hole pair 

is equivalent to a dynamically screened electron-electron scattering event because of the detailed balance between 
plasmons and electrons. Since ( p7| ) conserves total electron momentum (in the absence of Umklapp processes), its 
effect on the electrical current in materials with spherical energy surfaces is nonvanishing only because it rearranges 
the momentum distribution, which influences other scattering processes. The most important difference between 
using (^) and a previous treatment of the effects of electron-plasmon scattering on mobility is that the earlier 
work assumed an equilibrium (Bose-Einstein) distribution function for the plasmons. In the latter case, there is a net 
momentum loss from the nonequilibrium electron system. 



B. Electron-phonon collision term 

If the phonons interacted only with electrons by non-Umklapp processes, they would reach a local detailed balance 
with the electrons just as the plasmons do. But, in fact, the phonons interact with other phonons through anharmonic 
lattice forces. The excess momentum given to the LO phonons by the nonequilibrium electron system is dissipated 
primarily through interactions with acoustic phonons, which can lose momentum to the environment through umklapp 
processes. Generally, mobility calculations have employed the "Bloch assumption" that the phonon system may be 
treated as if it were in thermal equilibrium | p6[ | , valid when anharmonic processes are much faster than electron-phonon 
scattering. 

Adopting the model of equilibrium phonons, the electron-phonon component of W"^'^^ in (|3^) gives the collision 
term, 



-{/ (k)}^" 



coll 



±y±q 



[N\uJk,k-<i) + 1] lrnD{q, L^fe,fe_,)/(k) [1 - /(k - q)] 



Un"^ ^ I €u{q,UJk,k~q) 

-N°iLOk,k-g)lmDiq,LOk^k-g)fO^ - q) [1 - fM | • (48) 

The electron-phonon matrix element is screened by only the electronic part of the total dynamic dielectric function. 
The long-wavelength, static limit of the RPA e{q,uj) is the temperature-dependent Thomas-Fermi dielectric function 
with the high-frequency dielectric constant eoo- Therefore, equation ( p8[ ) in this limit gives a somewhat different result 
than treatments of electron-LO phonon scattering using Thomas-Fermi screening with the static dielectric constant 
eo |E2|. For the nondegenerate case, the long wavelength limit of (El) is in agreement with the original treatment of 



screening in polar semiconductors by Ehrenreich |27 

The ion-ion Coulomb interactions included in Hg (equation (^) make a self-energy contribution to the phonon 
Green's function |10|. The real part of this contribution shifts the LO phonon energy from hujTO to Klolo in the 
denominator of the noninteracting Green's function D^^"> — —2ujto / ^{i^lo ~ ^^)- Damping due to anharmonic 
interactions and disorder is neglected here. 

The interacting phonon Green's fimction 

- 2wTo/?i 

^'^'"'^ " 1 - M2Z?(0)x(g,a;) " c.^ - u^l^ - 2wToM^x{q,i^)/n ^ ' 

contains the self-energy term = M^D^^^'x due to the polarization of the electron gas. The real part of H'^^ 
renormalizes the LO phonon frequency, giving rise to two hybrid modes of mixed electron-phonon character. The 
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coupled modes have been treated previously with the finite-temperature plasmon-pole approximation for x [0 • 
In this approximation, the electronic excitation spectrum which couples to the phonons is represented by a single 
mode with energy varying between htUp = ^ Aime'^ /ni*eoo at small q and h?q^ /2m* at large q. The plasmon coupling 
strength and energy are determined by requiring fulfillment of the f-sum rule and zero-frequency Kramers-Kronig 
relation, so that the plasmon-pole approximation for x is 



X''{1.-)-^-^, (50) 



i^^'p[l-e-\q,Or- (51) 

Using the plasmon-pole susceptibility in the phonon self-energy H'^p gives a phonon Green's function with poles 
at the two frequencies w+ and uj^ , 

-I = ^ {-io + ± Mo - ^l? + H^^lo - -lo)] ■ (53) 

In the plasmon-pole model for the coupled modes, the phonon spectral function — 7rIm[Z3(g, w)] appearing in the 
electron-phonon collision term (^) has delta function peaks at ujj^ and 

Im[L>(g,w)] = Im[L'+(g,w) -hD-(q,tj)] (54) 

Im[I?±(g,a;)] = T^^y^^ + ^±) ~ S{lu - c.±)] . (55) 

nuj±[uj:^ — uj±) 

At low densities, the low-energy hybrid mode has frequency close to the uncoupled plasmon frequency ujp, but 
has nonzero weight in the phonon spectral function determined by the factor (cj^ — Ldp)/{uj'^ — uj'^). Even though 
this factor is small at low densities, the low-energy mode has an exponentially larger thermal occupation iV°(u;_) 
at low temperatures than the high-energy hybrid mode or uncoupled phonon mode do. Therefore, for low densities 
and temperatures, the presence of N*^ in ( [48|) suggests that the mobility calculated by including scattering from the 
coupled modes should be lower than the mobility calculated by using the uncoupled mode approximation. In addition, 
"antiscreening" effects due to the dynamic nature of the system response can generate screened potentials that are 
actually greater than in the unscreened case . Reference presents numerical results for n-type GaAs which show 
that including dynamic screening, mode coupling, and electron-electron scattering does significantly lower calculated 
mobilities in situations where ionized impurity scattering does not dominate, as in modulation-doped structures. 



IV. CONCLUSION 



The total dielectric function eT{q,i-o) lends itself to a simple and general method for deriving the inelastic collision 
term in the electron Boltzmann equation for scattering from a coupled electron-phonon system. When the Born 
approximation is valid, the inelastic differential scattering rate VF'"^' can be expressed in terms of the nonequilibrium 
total dielectric function ct, which includes screening by both electrons and phonons. In the RPA, W™"^^ separates into 
a purely electron-electron interaction plus an electron-phonon term describing interactions with the hybrid modes. 
These interactions are dynamically screened by only the electrons. Reference investigates the effects of the phonon 
self-energy, dynamical screening, and electron-electron scattering on the mobility in n-type GaAs by means of a 
numerical solution of the linearized Boltzmann equation. 

Because the total dielectric function approach to the Boltzmann equation treats the electrons and ions in a unified 
manner, it can be applied quite generally to a wide variety of coupled systems. More accurate and complete treatments 
1^,^ of the lattice part of e^^ could be useful for more complex electron-phonon systems. Also, it would be very 
interesting to try to include electron-impurity scattering through the 6^(5,0;) for an imperfect crystal with local 
phonon modes. Another logical direction to pursue would start from a dielectric function which goes beyond the 
RPA to include electron exchange effects in order to derive the exchange scattering term in the Boltzmann equation. 
However, a systematic approach to improving on the RPA is required, since exchange corrections can be largely 
canceled by correlation terms. 
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